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FOREWORD 


This  report,  "The  USAF  Stability  and  Control  Digital  Datcom,"  describes 
the  computer  program  that  calculates  static  stability,  high  lift  and  control  , 
and  dynamic  derivative  characteristics  using  the  methods  contained  in  Sec¬ 
tions  4  through  7  of  the  USAF  Stability  and  Control  Datcom  (revised  April 
1976).  The  report  consists  of  the  following  three  volumes: 
o  Volume  I,  Users  Manual 

o  Volume  II,  Implementation  of  Datcom  Methods 
o  Volume  III,  Plot  Module 

A  complete  listing  of  the  program  is  provided  as  a  microfiche  supplement. 

This  work  was  performed  by  the  McDonnell  Douglas  Astronautics  Company, 
Box  516,  St.  Louis,  MO  63166,  under  contract  number  F33615-77-C-3073  with  the 
United  States  Air  Force  Systems  Command,  Wright-Patterson  Air  Force  Base,  OH. 
The  subject  contract  was  initiated  under  Air  Force.  Flight  Dynamics  Laboratory 
Project  8219,  Task  82190115  on  15  August  1977  and  was  effectively  concluded 
in  November  1978.  This  report  supersedes  AFFDL  TR-73-23  produced  under 
contract  F33615-72-C-1067 ,  which  automated  Sections  4  and  5  of  the  USAF  Sta¬ 
bility  and  Control  Datcom;  AFFDL  TR-74-68  produced  under  contract  F33615-73- 
C-3058  which  extended  the  program  to  include  Datcom  Sections  6  and  7  and  a 
trim  option;  and  AFFDL-TR-76-45  that  Incorporated  Datcom  revisions  and  user 
oriented  options  under  contract  F33615-75-C-3043.  The  recent  activity  gener¬ 
ated  a  plot  module,  updated  methods  to  incorporate  the  1976  Datcom  revisions, 
and  provide  additional  user  oriented  features.  These  contracts,  in  total, 
reflect  a  systematic  approach  to  Datcom  automation  which  commenced  in  Feb¬ 
ruary  1972.  Mr.  J.  E.  Jenkins,  AFFDL  FGC,  was  the  Air  Force  Project  Engineer 
for  the  previous  three  contracts  and  Mr.  B.  F.  Niehaus  acted  in  this  capa¬ 
city  for  the  current  contract.  The  authors  wish  to  thank  Mr.  Niehaus  for  his 
assistance,  particularly  in  the  areas  of  computer  program  formulation,  imple¬ 
mentation,  and  verification.  A  list  of  the  Digital  Datcom  Principal  Investi¬ 
gators  and  individuals  who  made  significant  contributions  to  the  development 
of  this  program  is  provided  on  the  following  page. 

Requests  for  copies  of  the  computer  program  should  be  directed  to  the 
Air  Force  Flight  Dynamics  Laboratory  (FGC).  Copies  of  this  report  can  be 
obtained  from  the  National  Technical  Information  Service  (NTIS) . 

This  report  was  submitted  in  April  1979. 
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SECTION  1 


INTRODUCTION 

Digital  Datcom  calculates  static  stability,  high-lift  and  control 
device,  and  dynamic-derivative  characteristics  using  the  methods  contained  in 
Sections  4  through  7  of  Datcom.  The  computer  program  also  offers  a  trim 
option  that  computes  control  deflections  and  aerodynamic  data  for  vehicle 
trim. 

Even  though  the  development  of  '  Digital  Datcom  was  pursued  with  the 
sole  objective  of  translating  the  Datcom  methods  into  an  efficient,  user- 
oriented  computer  program,  differences  between  Datcom  and  Digital  Datcom  do 
exist.  Such  is  the  primary  subject  of  this  volume.  Implementation  of  Datcom 
Methods,  which  contains  the  program  formulation  for  those  methods  in  variance 
with  Datcom  methods.  Program  implementation  information  regardirg  system 
resources  necessary  to  make  the  program  operational  ate  presented  in  Sections 

5  and  6.  / 

/ 

/ 

Section  6  also  lists  each  of  the  routines  and  references  t^eir  appear¬ 
ance  in  the  program  listings  provided  as  a  microfiche  supplement  to  this 

/ 

volume. 

Users  should  refer  to  Datcom  for  the  validity  and  limitations  of  methods 
involved.  However,  potential  users  are  fore-warned  that  Datcom  drag  methods 
are  not  recommended  for  performance.  Where  more  than  one  Datcom  method 
exists,  the  summary  in  Table  1  indicates  which  method  or  methods  are  employed 
in  Digital  Datcom.  Tables  2,  3,  and  4  define  the  basic  output  data  in  each 
Mach  regime  and  shows  the  overlay  in  which  each  is  computed. 

The  computer  program  is  written  in  Fortran  IV  for  the  CDC  Cyber  175. 
Through  the  use  of  overlay  and  data  packing  techniques,  core  requirement  is 
67,000  octal  words  for  execution  with  th?  NOS  operating  system  using  the  FTN 
compiler.  Central  processor  time  for  a  case  executed  on  the  NOS  system 
depends  on  the  type  of  configuration,  number  of  flight  conditions,  and 
program  option  selected.  Usual  requirements  are  on  the  order  of  one  to  two 
seconds  per  Mach  number. 

Direct  all  program  inquires  to  AFFDL  FGC,  Wright-Patterson  Air  Force 
Base,  Ohio  45433.  Phone  (513)  255-4315. 
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Table  1  SUMMARY  OF  DIGITAL  DATCOM  METHODS 

AERODYNAMICl  DATCOM  MACH  METHOD!  T 

PARAMETER  CONFIGURATION  SECTION  REGIME  NUMBER  OVERLAY  SUBROUTINE  REMARKS 
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Table  1  SUMMARY  OF  DIGITAL  DATCOM  METHODS 


NDM-NO  DATCOM  METHOD  NP-NOT  PROGRAMMED 


Table  1  SUMMARY  OF  DIGITAL  DATCOM  METHODS 


NDM-NO  DATCOM  METHOD  HP-NOT  PROGRAMMED 


TABLE  3  OVERLAYS  DEFINING  EACH  OF  THE  BASIC  TRANSONIC  OUTPUT  PARAMETERS 
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12  •  SECOND  LEVEL  METHODS.  OVERLAY  35 


POWES  INCREMENTS 


SECTION  2 


PROGRAM  ORGANIZATION 

The  Digital  Datcom  program  consists  of  a  MAIN  program,  EXECUTIVE, 
subroutines,  METHOD  subroutines  and  UTILITY  subroutines.  The  organiza¬ 
tion  and  interfaces  between  these  program  components  are  shown  ir.  Figure  1. 
The  MAIN  program  performs  executive  functions  that  control  and  direct  all 
computations;  the  EXECUTIVE  subroutines  perform  noncomputational  tasks,  which 
include  input  data  manipulation  and  selection  of  output  formats;  UTILITY 
subroutines  perform  standard  mathematical  computations;  and  METHOD  sub¬ 
routines  implement  the  Datcom  stability  methods.  .. 
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EXECUTIVE  Subroutines 


FIGURE  1  OVERLAY  PROGRAM  STRUCTURE 


.SECTION  3 


EQUATIONS  FOR  GEOMETRIC  PARAMETERS 


One  of  the  main  features  of  Che  Digital  Datcom  program  is  that  a 
minimum  of  input  data  are  required.  Minimal  inputs  require  the  program 
to  calculate  basic  geometric  parameters  required  by  the  Datcom  methods. 
Equations  for  pertinent  geometric  parameters  are  defined  in  this  section. 
3.  1  PLANFORM  PARAMETERS 

The  nomenclature  used  in  the  equations  for  calculating  theoreti¬ 
cal  and  exposed  planform  areas,  taper  ratios  and  aspect  ratios  are  shown  in 
Figure  2.  Equations  for  these  parameters  are  presented  below  for  a  double 
delta  or  cranked  planform.  Straight-tapered  planform  parameters  are  obtained 
by  setting  b*  /2  ■  0.0,  C.  ■  C  ,  A*  »  1.0  in  the  following  equations: 

O  D  C  O 


b^/2  -  b/2  -  b*/2 
b£/2  -  b*/2  -  b*/2 
*1*  <hb/2)/(hb/2) 

XI  *  VCr 

cr  *  Cr[>l  *  (l  -  V 

X‘l  ' 

lo  '  Ct/Cb 

>*.»*»* 

w  I  o 

X  »  C  /C  ■ 
v  t  t 

»‘t  •  K  *  v  \/2 


SI  -  «V  *  V  V2 
<  -  <S  +  ct>  \n 
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SW  -  <Cr  +  V  V2  +  So 

A*  -  4(^/2)  2/sJ 

A*  -  4(b^2)2/S; 

Aj  -  4(b*/2)2/S* 

Ay  -  4Cb/2)2/sw 

Datcom  methods  use  correlations  that  are  based  on  wing  sweep  angles 

measured  at  various  chordlines.  The  nomenclature  used  to  calculate  sweep 

angles  is  presented  in  Figure  3.  Sweep  angle  equations  are  presented  below 

for  a  double  delta  or  cranked  wing.  To  obtain  straight  taper  wing  sweep 

angles  set  C  and  A  *  0  in  tne  following  equations: 

o 

Ct  -  4(1  -  X*^/£A*IC1  +  X*^] 

C0  "  4<2  ’  xV/tA*o(1+ XVJ 

Anj  ■  tan  2[Cj(n-n)  +  tanAm^] 

An  •  tan  2[C  (o-n)  +  t&nAa  ] 
o  o  o 

*  —1  *  *  *, 

(Ajj)ef:f  -  cos  ](Sj  cos  Anj  +  SQ  cos  AnQ)/S  ] 


The  nomenclature  used  to  calculate  the  exposed  mean  aerodynamic  chord 
(MAC)  for  a  double  delta  or  cranked  wing  is  shown  in  Figure  4.  The  param¬ 
eters  necessary  to  define  the  lateral  and  longitudinal  location  of  the 
exposed  MAC  are  included.  Equations  to  calculate  and  locate  the  MAC  are 
presented  below;  To  obtain  values  for  a  straight-tapered  wing  set  C*q  ■  0, 
Y*  *  6,  S*  *0  in  the  equations  below: 


— *  *  *  *9.  .  „ 

Cj  -  2Cr(l  +  Xj  +  Xj  )/3(l  +  Xj) 

.  *  *2..  .  * 

C  -  2C.  (1  +  \  +  \  )/Hl  +  X  > 

O  D  O  O  O 
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•m.  M  •VJ\ 


REQUIRED  INPUTS 
ALL PLANFORMS 

/''''  REQUIRED  INPUTS 

V  /  DOUBLE  DELTA 
'***  -<*  . 

AND 

CRANKED PLANFORMS 


m»  PERCENTAGE  CHORD 
AT  WHICH  SWEEP 
ANGLE  IS  DEFINED 

n  -  ANY  CHORD  LOCATION 
EXPRESSED  IN 
PERCENTAGE  CHORD 


FIGURE  3  SWEEP  ANGLE  NOMENCLATURE 


FIGURE  5  THEORETICAL  OR  REFERENCE  MEAN  AERODYNAMIC  CHORD 
NOMENCLATURE 


GLOVE  COMPONENT 


i 


T" . 1  ■  k 


BASIC  WING  COMPONENT 

TRAILING  EDGE  EXTENSION  COMPONENT 


*  -  c,- 
V  -  Ct  +  (bo/2)' 

b  12 
o 

<y  -  (s*  +  ct>  (bo/2)v 


(Ao)'  =  4[(b*/2)2]7(S*)' 


(X  )’  =  C  JC  ' 
o  t  fc 


Supersonic  nonstraight  wing  analyses  require  the  wing  to  be  synthesized 
from  basic  wing,  glove,  and  trailing  edge  extension  components  as  shown  on 
Figure  7.  When  the  leading  edge  outboard  sweep  angle  is  greater  than  the 
leading  edge  inboard  sweep  angle,  an  additional  geometric  parameter,  S^,  is 
required  and  is  shown  in  Figure  8.  Equations  for  calculating  geometric 
parameters  for  the  various  wing  components  as  required  by  the  stability 
methods  are  presented  below: 


All  Planfortns 

*  * 

<Cr*>bw  ■  S  +  <  r  -  TT-  3 11*"  Ve  -  *TE  1 

o  o 

*  KC*)  +  C  ]  b* 
basic  S,  -  lv  Vbw  V 

wing  &W  2 

component  *  ^*2 

\w  S* 

bw 

<C?>bw 


(C*>  -  (Tan  ATT?  )  ) 


glove 

component 


s*-  (C.)g  (f  -£) 

j2 

A*  -  l2  2  J 

8  ~n - 

v° 


45 


trailing 

edge 

extension 

span 


2  (—  -  k*  ) 

'2  ~7  J 


c*  fb* 

S%  *  {-— 


F  I  (>»->  -LE  ) 


S1  ’  5bw 


Geotaetric  parameters  required  for  horizontal  and  vertical  tail  anal>ses 
are  identical  to  those  for  wings.  Tail  parameters  can  be  calculated  by 
substituting  tail  geometry  for  wing  geometry  in  the  wing  Equations.  Vertical 
tail  lateral  stability  calculations  require  additional  geometry  parameters  as 
shown  in  Figures  9a  and  9b.,  Equations  are  listed  below: 


Straight  Tapered  Vertical  Tail 


C  -  C  -  (C  -  C)(ZJ/(b.  /2) 

v  r  r  i  u  v 

X  ■  *u  +  <V  -  \  -  ZH  <T“  \e  > 

V  I 

Non-Straight  Vertical  Tail 

n  zH  >  bv  b; 

H  ~2  "  T 

b  b* 


b*  b 
o  v 


+<V  Xv  -  <2 - 2>  XE,*  '  (ZH  +  r  '  T  >  TAa  *u, 


Cv-Cc  +  <S,-Ct>(r- V/(r> 

b  b* 

If  ZH  <  _v  _o 

“  2  “  2 

*  ■  \  *  b.  '  \  '  zh  'Veo> 

Cy  ■  Cr  -  (Cr  -  V(ZE>/  (\  .V,- 


For  a  horizontal  lifting  surface,  an  equivalent  dihedral  is  defined  as 


follows: 


r,  ,bf  *  .  „  ,b*  . 

eq  ’  ~2  0  r  T o 


6  x.4 


3.2  BODY  PARAMETERS 


Longitudinal  stability  analyses  for  bodies  in  the  supersonic  and  hyper¬ 
sonic  speed  regimes  require  the  body  to  be  synthesized  in  nose,  afterbody, 
and  tail  segment  components  as  defined  in  Figure  11.  Geometry  parameters  for 
the  various  body  segments  analyses  are  defined  below: 


lB  +  *A 

lBT  iB  “  lB 


S 

P 


rx  (dl,) 


T' « «■» 

^  o 

Sp 

f0ib  SX 


Body  planform  area 

Body  base  area 

Distance  from  nose  of  body  to  centroid  of 
planform  area 

Volume  of  body 


If  d 2  >  d^,  calculate  flare  angle  0^  ■  TAN 

If  calculate  boattail  angle  ■  TAN  *  |*^^1  2^ 

L  iBT 


3.3  GENERAL  SYNTHESIS  PARAMETERS 

Synthesizing  and  interference  nomenclature  for  longitudinal  and  lateral 
stability  calculations  are  defined  in  Figure  12.  The  geometric  parameters  are 
presented  in  equation  format  below; 

-  (b/2  -  b*/2)  TAN  AOj  COS  (o^ 

“c,  •  *c„  -  *  “V 

(X4C)W  ■  (Xac/C£)w  CJ;  Where  (X4<,/C*)  is  calculated  in  wing  pitching 


moment  subroutine 


POSSIBLE  SUPERSONIC  AND  HYPERSONIC  BODY  CONFIGURATIONS 


in 


in 

Ik 

iBT*° 


S 

*1**2 


NOTES: 

NOSE  AND  TAIL  SEGMENTS  MAY  BE  CONICAL 
(AS  SHOWN)  OR  OGIVAL 

DIAMETERS  dH,di,  AND  d*  ARE  COMPUTED 
FROM  LINEAR  INTERPOLATION  OF 
INPUTS  t-  VSR 


FIGURE  11  SUPERSONIC  AND  HYPERSONIC  BODY  GEOMETRY 


(AX  )  -  AX  -  (X  )  COS  (a  ) 

ac  w  eg  ic  v  l  w 

AXjj  -  (b/2  -  b*/2)H  TAK  AOj  COS  (ct^ 

H 

<4Xcg>H  •  Xcg  '  (XH  *  4XH> 

2*  ■  2h  -  TAS  (at)H 

(Xac>H  '  <Xac/C?>H  C‘v 

<Zac>H  *  ZH  -  <Xac>H  518  <»1>H  ’  Z=g 

4<Xac>H  ■  <1Xcg>a  -  <Xac>H  C0S  (ai>H 

<XC/4>K  '  *H  '  (Xr>H  C0S  <“l>H 

Z'  •  -Z  +  (C  /A)  SIN  a. 
w  w  r  i 

*f  ■  x.  +  4X»  ♦  &  >  148  4le  +  T“  \.e,  ♦  ^ 

2  0  2  12 

i  -  A  -  X  +  (X  )  +  (C  ) 

p  V  eg  r'w  rv 

4 

Zp  -  Z  +  (YJ 
K  eg  R  v 


3.4  DOWNWASH  PARAMETERS 

Oownwash  geometric  nomenclature  is  defined  in  Figure  13.  The  equations 
presented  below  are  used  primarily  in  the  subsonic  speed  regime: 

ZH  “  ZH  "  \  SIN  (ai)H  “  Zw  ♦  c  SIN  <ai) 

H  rw  W 

,  ni  *  h  +  Xru  C0S  <oi>H  -  (Xw  *  Cr  C0S  {ai)J 
H  w 

A4  -  Z’  TAN  (ai)w 

Lt  -  lH  -^Lh 

AhHi  -  z’/COS  (ai)w 
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reference  < 
plane 


*h2  -  Lt  SIN  (oi)w 

hH  *4\  +4hH2 

l 0  -  L  COS  (ai) 

2  T  w 

Y  -  ARCTAN  (hn/i2) 

*3  *  <Cr>„  -  <Vw 
If  beff/2  <  (b/2  -  b*/2) 

C 


o  w 


Ceff  ’  Crw  "  b^2  -**f2  (beff/2) 


w 


ha  -  (beff/2)  TAN  A°i  +  ct  « 

eff 


(Eeff  -  C.  ) 


eff 


w 


If  beff/2  9  (b/2  ~  bS/2) 


eff 


c  -  c 

S  -  bi/2  ‘  lbeff/2  -  <b'2  - 


Eef f  "  (b/2  '  bJ/]>w  TAN  +  tbeff/2  -  (b/2  -  b*/2)w]  TAN  Aoq  +  C  /4 

eff 


*eff  "  *2  "  (Eeff  "  Cr  J 

w 


3. 5  POWER  EFFECTS  PARAMETERS 

Geometric  parameters  required  to  calculate  propeller  and  jet  power 
effects  are  defined  in  Figures  14  through  18.  Power  effects  are  only  cal¬ 
culated  for  longitudinal  stability  results  In  the  subsonic  speed  regime. 
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POWER  ON 


PROPELLER 


MAC  ^-HORIZONTAL 
4  J  TAIL  MAC 


f  -ZJ  /-WING  MAC 

+X 

\ 

MAC 

+% 

« 

Xp 

_xp 

■ 

1* 

=  ZV-XT  SIN  a: 

"  'w  'w 

i 

“p 

=  <*SCH  +ait+*u"*p 

z* 

=  ZT  +Xp  TAN  a  p 

1 

\ 

=  Zh-ZT+[{Xt,+XrhCOS 

Ib-*P»TM«|T| 

4 

|<-’V0S%| 

\ff 

“W'hTAIMa'p- 

‘  POWER* 

FIGURE  14  DEFINITION  SKETCH  FOR  PROPELLER  POW 

:R  EFFECT  CALCULATIONS 

SINGLE  ENGINE 


FIGURE  15  GEOMETRY  FOR  DETERMINING  IMMERSED  WING  PARAMETERS 


SINGLE  ENGINE 


FIGURE 16  GE0METRY  F0R  DETE™'-^  immersed  W,NG  PARAMETERS  (CONT’D) 
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*K3  *&*■  >_^t®»to«jvV  (Wfclii 


2 


IFb*0/2=  0.0 
IF  yp  =  b/2  -  b*0/2 
IF  yp  >  b/2  -  b*0/2 

S*j=  2[2(bj/2)]  Crp 

[2(bi/2)]  2 
A*i=  0.50  S*j 

A*  =  1-0 


FIGURE  17  GEOMETRY  FOR  DETERMINING  IMMERSED  WING  PARAMETERS  (CONCLUDED) 


VCf'  b/2 

(Cr  -  C„)  yt, 

1  ~  b/2-b*^ 


VC*~ 


(C„  -  Ct)  (yp  -  b/2  +  b'j/2) 


bV  . 

IMMERSED  AREA  FOR  TWO  ENGINES 
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.  I 


*  ’•  / 

t 


cos«it 

2  j  *  (XH+  ( Y  COS«j()  -  Xe)  SIN ♦  (Zir  ZT)  COS^ 

X'j »  4.6  Rj 

*i.**j**'. 

FIGURE  18  DEFINITION  SKETCH  FOR  JET  POWER  CALCULATIONS 
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3.6  GROUND  EFFECTS  PARAMETERS 

Ground  effects  are  only  calculated  for  longitudinal  stability  results  in 
the  subsonic  speed  regime.  Lifting  surface  heights  that  are  required  by  the 
Datcon  ground  effect  analyses  are  defined  In  figure  .9  and  are  presented  in 

equation  format  as  follows: 

Equations  for  Calculating  hc./5b/2 

IF  rt  *  0  ANO  ib  ?)f  =  0.25  (b/2)  h0.75b  7  *  H.ttCi  * UK 

IF F j *  0  ANO  (b  2) r ^  >  0.25  (b 2)  "C^C,3  TANro[(b7Y#  “Slb/aj+AXTAMI.ft 

IF  r :  ^  0  AND  (b  7)p  *  0.25  (b  7)  N)  75b  7  3  **0 ,75Ct  *  0  75  (b 71 TAH  ri*  AX  ™  <*'V| 

#  I  I 
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Equations  foe  Calculating  h 
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«  I 
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Equations  for  Calculating  H 
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Equations  for  Calculating 


(hcf<)H=HG  +  zH-^HTAN<VH 


IF  Tj  =  0  AND  fyf»H  =  (b  2>H  -  <b  ?y 
"  i  oHl 

IFriH=OANDfyt)H  |(b^H-(b?VoJ 
IFf-jf  OANO(7,'h  •  [(b,2)H  "  ,b  7V 


IFT^  OANO(7()h  |(b.^H-(b/2)r 
I  °H 


T,Br" 


Ground  effect  methods  require  calculation  of  a  planform  parameter 
AX,  in  addition  to  the  previously  defined  ground  heights.  This  param 
eter  is  shown  in  Figure  20. 


Straight  Tapered  Wing 

AX  =  0.75  C,  -0.75  (fa,®  TAM  A*a 

Cranked  or  Double  Delta  Wing 

IF  bo72  “  °-25  (b/2V  AX  =  0,75  Cf  “  °-75(b/?)TAN  A25j 

IF  b*0/2  >  0.25  AX  =  0.75  Cr  -  TAN  ^ | -  0.25(b/2)J  -  TAN  A25i |(b/2)  -  b 


Straight  Tapered  Wing  Cranked  or  Double  Delta  Wing 


FIGURE  20  GROUND  EFFECTS  PLANFORM  PARAMETER  AX 


SECTION  4 


INCORPORATION  OF  METHODS 

This  section  summarizes  those  methods  which  were  incorporated  into 
Digital  Datcom  but  were  not  defined  in  the  Datcom  Handbook  or  involve  method 
interpretation.  Though  some  of  the  methods  included  are  not,  in,  general, 
standard  Datcom  methods,  they  permit  greater  flexibility  in  using  the  pro¬ 
gram,  and  provide  output  for  some  parameters  which  can  be  closely  approxima¬ 
ted  or  are  difficult  to  obtain  experimentally.  All  of  the  methods  presented 
in  this  section  are  referenced  to  Table  1  of  Section  1  and  the  Datcom. 
Methods,  or  procedures,  not  outlined  in  this  section  follow  the  Datcom  method 
and  users  should  consult  the  Datcom  for  method  limitations  and  formulation. 
4.1  AIRFOIL  SECTION  AERODYNAMICS 

This  section  describes  a  procedure  that  can  be  used  to  obtain  the 
geometric  and  aerodynamic  section  characteristics  of  virtually  any  user 
defined  airfoil  section.  Its  incorporation  into  Digital  Datcom  frees  the 
user  from  the  labor  of  calculating  those  section  parameters  that  were  re¬ 
quired  Inputs,  yet  allow  him  the  flexibility  to  alter  those  parameters  for 
which  he  has  data. 

The  Airfoil  Section  Module, will  accept  the  following  user  inputs: 
o  The  airfoil  section  designation 
o  Section  upper  and  lower  cartesian  coordinates 
o  Section  mean  line  and  thickness  distribution 
By  these  three  methods, ,  many  airfoil  sections  can  be  described  and  the  sec¬ 
tion  characteristics  calculated. 

Since  the  Airfoil  Section  Module  (ASM)  uses  the  Mach  and  Reynolds  number 
inputs,  they  must  be  defined  in  namelist  FLTC0N  using  MACH  and  RNNUB.  How¬ 
ever,  the  ASM  uses  the  unit  Reynolds  number  and  by  implication  treats  a 
section  one  foot  (or  meter)  in  length. 

This  module  brings  together  the  outstanding  features  of  two  separate 
studies.  Kinsey  md  Bowers  (AFFDL-TR-71-87)  have  written  a  program  that 
calculates  the  airfoil  coordinates  of  select  NACA  designations,  then  uses  the 
Weber  technique  to  calculate  the  section  aerodynamic  characteristics. 
Nleldling  of  McDonnell  Aircraft  has  written  a  similar  program  using  the  Weber 
method,  then  incorporates  additional  methods  to  refine  the  theoretical 


TABLE  5  AIRFOIL  SECTION  MODULE  ROUTINE  DESCRIPTION 


PROGRAM/SUBROUTINE  PURPOSE 


M500?2  (OVERLAY  50.0) 
INIZ 
SECI 
SECO 
CSLOPE 
XYCORO 
OELY 

AIRFOL  (OVERLAY  (50,1)) 
DECODE 
C00R04 
C00RD4M 
COORDS 
COCRD5M 
COOROI 
COORD6 
COR DSP  , 

SLEQ 

THEORY  (OVERLAY  (50.2)) 
IDEAL 
SLOPE 
ASMINT 


MODULE  EXECUTIVE  PROGRAM 
INITIALIZE  IOM 
READ  USER  INPUTS 
TRANSFER  MODULE  OUTPUTS 

CALCULATE  VARIABLE  SLOPE  FOR  SUPERSONIC  AIRFOILS 
CALCULATE  AIRFOIL  SECTION  FROM  USER  INPUTS 
CALCULATE  DATCOM  PARAMETER  AY 

MAIN  PROGRAM  FOR  NACA  DESIGNATION  INPUTS 
RE  AO  USER  INPUT  NACA  DESIGNATION.  DECODE 
CALCULATE  4-DIGIT  NACA  AIRFOIL 
CALCULATE  4-DIGIT  (MODIFIED)  NACA  AIRFOIL 
CALCULATE  5-DIGIT  NACA  AIRFOIL 
CALCULATE  5-DIGIT  (MODIFIED)  NACA  AIRFOIL 
CALCULATE  1-SERIES  NACA  AIRFOIL 
CALCULATE  6— SERIES  NACA  AIRFOIL 
CALCULATE  SUPERSONIC  AIRFOIL  COORDINATES 
SIMULTANEOUS  LINEAR  EQUATION  SOLVER 

MAIN  PROGRAM  FOR  AIRFOIL  AERODYNAMICS 
CALCULATE  SECTION  IDEAL  AERODYNAMICS 
CALCULATE  LIFT  AND  MOMENT  SLOPES 
NON-LINEAR  INTERPOLATION  ROUTING 

CALCULATE  VARIABLE  CLMAX  FOR  SECTION 


MAXCL  (OVERLAY  (50,3)) 


predictions.  A  cross  of  the  two  procedures  (coordinates  of  NACA  airfoils  and 
viscous  correction  from  Kinsey  and  Bowers,  and  the  aerodynamic  methods  of 
Nieldling)  yields  a  program  that  generates  fairly  accurate  results. 

The  module  is  incorporated  into  Digital  Datcom  as  Overlay  50,  and 
includes  three  secondary ; overlay  programs.  The  routines  use  the  10M  arrays 
for  data  storage  so  that  core  size  will  be  kept  to  a  minimum.  Table  5 
describes  each  of  the  22  module  routines  and  the  logic  flow  of  the  module  is 
presented  in  Figurees  21  through  24. 

4.1.1  Weber's  Method 

The  calculation  of  the  pressure  distribution  over  the  surface  of 
an  airfoil  in  an  incompressible  inviscid  flow  is  accomplished  by  use  of 
the  method  of  singularities.  Conformal  transformations  are  used  as  an 
intermediate  step  in  deriving  the  methods  for  determining  the  distributions 
of  singularities  from  which  the  velocity  distributions  are  calculated.  The 
routine  inputs  are  the  airfoil  coordinates  distributed  in  any  fashion,  the 
angle  of  attack,  and  the  Mach  number.  The  airfoil  shape  is  defined  by  curve 
fitting  the  input  coordinates  to  obtain  the  airfoil  geometry  at  thirty-two 


required  points,  i.e. : 


X  -  0.5  (COS  9v  +  1) 

0  *  vrr/32  for  o<v<32 


The  chord  line  is  obtained  by  joining  the  leading  and  trailing  edges 
of  the  airfoil,  where  the  leading  edge  is  defined  as  the  forward  most  point 
so  that  all  points  on  the  airfoil  surface  have  a  positive  x  coordinate. 

The  airfoil  is  placed  in  a  uniform  stream  Vq  at  an  angle  of  attack 
relative  to  the  chord  line.  The  velocity  Vq  is  resolved  into  components 
parallel  and  normal  to  the  chord  line. 

V  -  V  cos  a 
xo  o 

V  *  V  sin  i 

20  O 

Combining  the  results  for  the  parallel  and  normal  flows,  the  velocity 
distribution  equation  for  a  symmetrical  airfoil  at  angle  of  attack  is 

V  * 

o  (  r  1  r  dz  <*x' 

V(XfZ)  -  |  cos  a  |1  +  —  /  — ,  —■ , 

1  +  (dz/dx)  J 


l"“  [‘  +  ?/[ 


1  -  (1 


-2x')J  y 
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OVERLAY 

(50,0) 


FIGURE  22  AIRFOIL  SECTION  MODULE  -  NACA  DESIGNATION  ROUTINE 
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CALCULATE  IDEAL  AERODYNAMIC 
PARAMETERS 

ai>  ao*  aOj_<  ^Lj* 


CALCULATE  C»  @M  =  0 
Ja 


MACH  LOOP,  CALCULATE 

C  i  AT  EACH  MACH  NUMBER  INPUT 
Ja 


FIGURE  23  AIRFOIL  SECTION  MODULE  -  SECTION  AERODYNAMICS  ROUTINE 


.  '  . /  . . . 


In  the  Weber  Method  certain  combinations  of  the  above  terms  have  been  rede¬ 


fined  as  follows: 


SU)(x) 


d  z  d  x ' 
dx*  x  -  x' 


(Function  for  Source  Distri¬ 
bution  in  Parallel  Flow) 


S(2)(x) 


s(3)<*>  -i  /'fe - ^-t)^ 

\  1  -  (1  -  2x’)  ' 


(Slope  of  Thickness 
Distribution) 

•  (Function  for  Vortex  Dietri- 
bution  in  Normal  Flow,  to 
Account  for  a) 


These  functions  are  approximated  by  sums  and  products  of  the  airfoil  ordi¬ 
nates  and  certain  coefficients  which  are  independent  of  the  section  shape 


N-l 

S(1)(x  )  -  ^ 


s  (1)  ? 

vt-  rv 


S(2)(x) 


c  (2)  , 

*Vl>  V 


s<3)(x  )  =  2  s 


s  (3)  z  +  S' 

UJ/  V  Ni> 


V21 


The  effects  of  camber  on  the  resulting  velocity  distribution  are 

obtained  by  assuming  the  camber  to  be  small  compared  with  the  chord. 

This  results  in  the  camber  effect  being  accounted  for  in  the  parallel 

flow  V  *■  V  cos  a  only, 
xo  o 

The  Vortex  Distribution,  7(X),  on  the  chord  line  which  produces  a 
given  velocity  normal  to  the  chord  line  and  which  is  zero  at  the  trailing 
edge  is 


?(x„  ) 

2V._  ' 


s  (4)z 

VV  £ 


(Vortex  Distribution  due  to 
S  (xu  >  Camber) 


The  total  velocity  Vx(x,o)  on  the  chord  line  for  an  airfoil  with  camber  and 
incidence  is 

V  (x,0)  =  V  cos  a  \  1+  S  ^  (x) 
x  o  L 


with  the  +  sign  being  for  the  upper  surface  and  the  -  sign  for  the  lower 
surface. 


The  resulting  velocity  distribution  at  the  airfoil  surface  is  computed 
using 


(Slope  of  Camber  Line) 


where 


which  is  the  complete  expression  for  an  arbitrary  airfoil  at  angle  of  attack 
in  an  ideal  flow.  The  S^/j)(X)  and  (X)  terms  are  computed  by  approxima¬ 
tion.  The  pressure  coefficient  is  obtained  by 


The  term  1  +  S  (x 


(!),,,  .  AU) 

(x)  +  S  (x)  acc 


mints  for  the  vorticles  being  put  into 


a  flow  with  velocity  V  (1  .+  S^'(x)  +  S^^(x))  instead  of  V  .  The  term 

o  o 

(1  +  S*1  (x))  accounts  tor  the  differences  in  the  vortex  distribution  between 

the  thick,  and  thin  wing.  The  term  1/  (1  +  [Stfc/(x)  +  S'  tx)]"]  is  the  cor¬ 
rection  between  velocities  on  the  chord  line  and  on  the  surface. 


A . I . 2  Compressibility  Correction  and  Integration 

The  effects  of  compressibility  are  accounted  for  in  Weber's  Method  by 
the  application  of  compressibility  factors  to  the  velocity  distribution 
contributions  due  to  thickness  and  camber,  respectively. 


U-yi  -  m;> 


» *«lh  -  <>  -  <■„ )) 

1  o  o  P . 


O  +  s(1>,“ 
I  .  (s(->)2 


The  velocity  distribution  in  compressible  fle«  is  then  given  by 

I  r  S<»  S<4’1  sinof  S<»1  !-■ 


^  1 


i  s  s  s  i  n  a  s  ! 

)*'os  a  l  - —  +  — *  +  — —  1  +  r —  V 

!  i-  ~  $  J  -  0  L  bJ" 


The  compressible  pressure  coefficient  from  the  compret si ble  form  of 
Bernoulli 's  equation  is 


:hl . -t  4)  ''I]’-’ 


The  airfoil  lift,  axial  force  and  pitching  moment  are  computed  iroo  the 
compressible  and  incompressible  solutions  in' the  following  manner 


bet  I  ,  r  ;v>  .  r  (><) 

*  P.  ’P,  ■ 


"  r  /  /,  - 


W  if  W 


where  w  . 


.  V  •  ty 


7*, 


Therefore  trapezoidal  rule 

/  r 

CL  (M)  =  —  Mx 

N  (  . 
t r 

i'os  a 

Similarly 

13-1 

n  y 

CA  (M)  =  —  .f- 


«  JL  k 

|/x 

sin  9  1  ■+  ^x 

s  in 

■ 

j 

S-J 

«.  y 

lx 

sin  9  |  .  / 

N  i 

L 

9  1 

0 

•  2 

i^=i 

2  ■  K) 

>M  . 

y  ‘x  sin  9 


l'°SON  *  =  l  I  2  li/ 


•j-i  r  .  (2) 

n  y  f-  (SKi*(x)  + 

M)  =  —  h\  L 


S(5)(x)>  -  C  (M)  (S(2)(x) 

P1 


CM  (M)  * 


+  1/2  r  (M)  y2p 


-  r  [/ 


x  ( y- .  2  r>  1 


4.  1.  J  Ideal  Parameters 

Tlie  ideal  parameters  are  obtained  from  thin,  airfoil  theory,  which 
in  effect  means  results  are  obtained  for  tne  meanline  characteristics  in  an 
incompressible  inviscid  flow.  The  ideal  angle  of  attack  *  is  obtained  from 


t-lx 

a ,  =  I  /  - - 

i  /  r  • 

J.  ir  x(t-x 


r  1  i/2 

[“'H 


However,  at  the  leading  and  trailing  edges  the  equation  is  undefined  and 
increments  in  the  vicinity  of  the  leading  and  trailing  edges  must  be  deter¬ 
mined,  in  addition  to  the  Integration  over  the  interior  portion  of  the 
chord. 


X  “ 

X  * 


■*  .  i ;  pi  z 


<{/ 

dx 


n  to 

. iil 


x  *  .(MSI 


x*d 
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4.1.4  Crest  Critical  Mach  Number 


The  crest  critical  Mach  number  is  precisely  defined  as  that  free  stream 
Mach  number  for  which  local  sonic  flow  is  first  reached  at  the  3irfoil 
surface  crest  on  the  assumption  of  shock  free  flow.  Its  significance  is 
founded  on  its  relation  to  the  drag  rise  Mach  number.  Various  empirical 
studies  have  been  aimed  at  finding  the  critical  pressure  ratio  at  the  crest 
which  corresponds  to  a  drag  rise  in  the  test  data.  Nitzberg  (NACA  KMA9G20) 
proposed  a  critical  pressure  ratio  for  drag  rise  of 

pCREST/pTOTAL  =  U.  52«3 

which  corresponds  to  a  crest  Mach  number  of  M  »  1.0.  Sinnot  (RAS  TOM-6407) 
proposed  the  ratio 

pCRKST/pTOTAL  a  0.515 


which  corresponds  to  a  Mach  number  at  the  crest  of  M  *  1.02  and  which  corre¬ 
lates  better  with  drag-rise  data.  Sinnot's  value  is  used  in  the  Airfoil 
Section  Module,  thus  the  crest  critical  Mach  number  corresponds  to  a  local 
flow  at  Mach  1.02  at  the  crest  rather  than  sonic  conditions.  The  relation¬ 
ship  between  the  crest  pressure  and  crest  critical  Mach  number  is 


0.515(1  +  0.2  >I(2C)3'5  -  1 


CREST 


0.7  F  M 


CC 


where 


F  “  |*CC  +  3/2  U  -0rr>C. 


'CC' 


^CREST 


]- 


Mcc  -  CREST  CRITICAL  MACH 


C  -  INCOMPRESSIBLE  VALUE 

PCREST 


CC 
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Rewritten  so  that  Mcc  is  a  function  of  Cp(;REST,  the  relation  is  approxi¬ 
mated  by 

MCC-  [l.023  -  .9507  CpcR£ST  -.414  CpgREST  -  .1506  Cp^REST  -  .0212  C^rest] 

The  crest  location  for  each  angle  of  attack,  is  determined  by  comparing 
the  airfoil  surface  slope  for  each  x  location  to  tangent  u.  The  final 
location  is  obtained  by  interpolating  between  the  two  given  x  locations  whose 
airfoil  slopes  bracket  the  tangent  a  value.  The  CpcREsT  vaiue  is  obtained 
by  interpolation  of  the  Weber  incompressible  pressure  distribution  between 
the  two  x  values  surrounding  XcREST*  The  crest  critical  lift  coefficient 
is  obtained  using  the  Karman-Tsien  compressibility  rule  on  the  M  -  0  inte¬ 
grated  Weber  lift  coefficient. 


where,  CL(M)  -  CL  for  M  -  0. 

No  specific  boundary  layer  correction  is  used.  However,  the  Uatcom 
recommends  a  5Z  correction  factor  to  bring  the  results  in  line  with  experi¬ 
mental  data,  and  the  viscous  correction  of  section  lift  curve  slope  proposed 
by  Kinsey  and  Bowers  (Appendix  B,  Volume  I)  has  been  incorporated. 
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4 .2  TRANSONIC  WING  Ct.  FAIRING,  TRANSONIC  WING  FAIRING,  and  TRANSONIC 
WING  Cn  FAIRING 

Datcom  wing  methods  in  the  transonic  Mach  regime  calculate  aerodynamic 
parameters  only  at  specific  Mach  numbers.  Data  at  the  requested  Mach  number 
is  then  determined  by  interpolation.  This  approach  is  used  for  the  wing 
lift  curve  slope  (Ci  ),  wave  drag  (Cd-),  and  aerodynamic  center  (Xac). 

>  l  »V 

Nonlinear  fairings  for  each  of  these  parameters  are  discussed  in  the  follow¬ 
ing  paragraphs. 

4.2.1  Transonic  Fairings  of  Wing  C; 

Wing  lift  curve  slope,,  C^  ,  is  calculated  in  subroutine  TRS0NI,  overlay 
24.  The  same  methods  are  used  for  the  horizontal  tail  in  subroutine  TRS0NJ, 
also  in  overlay  24. 

Datcom  section  4. 1.3. 2  defines  the  methods  for  calculation  of  Cj  at 

a 

five  discrete  Mach  numbers  from  0.6  to  1.4.  Values  at  Mach  0.6  and  1.4  use 
the  subsonic  and  supersonic  methods,  respectively.  The  routine  used  to  fair 
this  curve  is  a  modified  version  of  subroutine  ASMINT  used  in  the  Airfoil 
Section  Module,  overlay  50.  To  ensure  a  smooth  continuous  interpolation,  a 
curve  is  constructed  by  fitting  the  points  by  a  left-hand  parabola  joined  to 
a  series  of  cubic  curves,  and  finally  a  right-hand  parabola.  This  technique 
yields  c  function  which  has  continuous  derivatives  everywhere.  The  slope 
of  the  curve  at  subsonic  Mach  numbers  is  obtained  by  differentiating  the 
equation  on  Datccm  page  4.1.3.2-49  with  respect  to  Mach  number.  At  Mach  1.4 
th»  slope  is  found  by  calculating  values  at  Mach  1.3,  1.4  and  1.5  and  assum¬ 
ing  a  curve  of  the  form: 

CL  -  A  +  B/tf  +  C/s- 

Subsonic  methods  are  used  to  Mach  0. 75,  or  0.1  less  than  the  force  break 
Mach  number  ( M f t> ) ,  whichever  is  smaller,  and  transonic  fairings  are  initiated 
at  that  point. 

Subroutines  TRANWG  and  TRANHT  are  used  to  calculate  at  Mach  1.3, 

1.4,  and  1.5  and  return  Cr  and  its  slope  at  Mach  1.4.  Subroutines  TRS0NI 

a 

and  TRS0NJ  calculate  using  the  subsonic  equation  if  the  Mach  number 

is  less  than  0.75  (or  Mf'j,  -  0.1),  calculate  the  slope  of  the  subsonic 
curve  at  Mach  0.75,  and  call  the  new  fairing  routine  if  the  Mach  number  if 
greater  than  0.75. 
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4.2.2  Transonic  Fairing  of  Wing  Cn^ 

The  wing  wave  drag,  C^,  is  calculated  in  subroutines  TRS0NI  and  TRS0NJ, 
overlay  24,  for  the  wing  and  horizontal  tail,  respectively.  The  method  is 
given  in  Datcora  section  4.1.5. 1. 

Digital  Datcom  performs  a  linear  interpolation  of  Datcom  Figure  4. 1.5.1- 
29  at  fifteen  discrete  Mach  numbers  to  determine  the  variation  of  .  Non¬ 
linear  interpolations  of  this  curve  are  performed  as  required  at  tiie  user 
defined  Mach  numbers  using  the  fairing  routine  developed  for  wing  .  Two 
additional  constraints  were  applied  to  perform  this  fairing. 

a.  If  the  linear  slope  to  the  left  or  right  of  a  given  point,  except 
the  end  points,  is  less  than  UNUSED,  (10  ^  on  CDC  computers),  the 
slope  at  that  point  is  set  to  zero. 

b.  Any  computed  value  less  than  zero  is  set  to  zero. 

Within  the  fairing  routine,  the  number  of  points  in  the  curve  is  used  to 

discriminate  between  a  fairing  of  Cn  and  C.  . 

w  L“ 

4.2.3  Transonic  Fairings  of  Wing  Aerodynamic  Center 

Aerodynamic  center,  X  ,  is  calculated  in  subroutines  TRANCM  and 

3  C 

TRHTCM,  overlay  25,  for  the  wing  and  horizontal  tail,  respectively. 

Datcom  section  4. 1.4. 2  defines  the  method  for  calculation  of  X  at 

ac 

six  discrete  Mach  numbers  from  0.6  to  1.4.  Values  at  0.6  and  1.4  are  deter¬ 
mined  using  the  subsonic  and  supersonic  methods,  respectively;  the  remaining 
four  points  are  obtained  from  Datcom  Figure  4.1.4.2-30  corresponding  to 
V  -  -2,  -1,  0  and  +1.  If  the  thickness  ratio  is  less  than  or  equal  to  IX, 
these  data  are  interpolated  for  the  aerodynamic  center.  If  the  thickness 
ratio  is  greater  than  IX,  the  curve  is  defined  Using  points  which  are  a  func¬ 
tion  of  the  force  break  Mach  number,  M^*'  An  increment  to  the  aerodynamic 
center  is  found  from  Datcom  Figure  4.1.4.2-33  and  applied  at  the  fifth  point 
(M^  +0.07)  arid  the  resulting  curve  is  then  interpolated  for  the  aerodynamic 
center.  The  following  table  summarizes  the  Interpolation  table: 
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Using  Six  Points 
t/c  <  7% 


Using  Eight  Points 
t/c  >  7Z 


M1 

0.60 

0.60 

M2 

M 

for  V  ■ 

-2 

(0.60fMfb)/2 

M3 

M 

for  V  “ 

-1 

Mfb 

M4 

M 

LOT  V  ” 

0 

Mfb  +  0.03 

*5 

M 

for  V  * 

+1 

Mfb  +  0.07 

“6 

1.40 

Mfb  +  0.14 

*7 

- 

M  for  V  -  +1 

M8 

1.4 

The  interpolation  routine  used  is  similar  to  the  routine  used  for  C.  and  Cn 

'  L  a  uu 

(Sections  4.2.1  and  4.2.2). 

4.3  TRANSONIC  WING  CL  ,  TRANSONIC  WING  Cp  ,  TRANSONIC  WING  •BODY-TAIL  Cp  x 
TRANSONIC  WING-BODY-TAIL  C„,  TRANSONIC  frlNG  CA  .and  TRANSONIC  WIn8- 
BODY  C„ 

This  section  describes  those  methods  used  to  compute  the  transonic  con¬ 
figuration  aerodynamics  using  Second  Level  Methods,  and  are  summarized  in 
Table  6.  Additionally,  the  partial  output  is  described. 

4.3.1  Transonic  Wing  Lift  Coefficient,  C. 

The  wing  lift  curve  versus  angle  of  attack  is  programmed  in  subroutine 
WINGCL.  The  method  described  in  Datcom  section  4. 1*3. 3  is  used  as  a  guide  to 
produce  trends  and  is  not  construed  to  be  an  exact  method  of  solution.  Since 
the  method  is  an  approximate  one,  the  following  procedure  was  employed  to 
produce  the  wing  lift  characteristics  applicable  to  thin,  low  aspect  ratio 
wings: 

1.  The  required  experimental  data  inputs  by  the  user  are  (zero 
lift  angle  of  attack)  and  aA  (the  angle  of  attack  where  the  lift 
becomes  nonlinear). 

2.  The  lift  variation  is  assumed  to  be  linear  up  to  aA,  and  nonlinear 

to  ar  (maximum  lift  angle  of  attack). 

L 

max 
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‘C°oWBT  IS  AVAILABLE  FR0M  THE  SECOND  LEVEL  ROUTINE  OF  OAT  COM,  SECTION  4.5.3.1, 
SUBROUTINE  WBTCD0. 
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3.  The  nonlinear  lift  region  is  modeled  by  a  mathematical  relationship 
that  satisfies  the  following  conditions: 


at  n  m  ci 


C,  (a*  -  a  )  at  a* 


at  a  *  a 


at  *  ft 


A  modified  polynomial  of  the  form 


A  +  B(X-X  )  +  C(X-X  ) 
o  o 


is  utilized  to  satisfy  each  of  the  boundary  conditions  and  yield  a  curve 

somewhat  parabolic  in  shape.  This  relationship  has  provided  excellent 

results  in  modeling  the  nonlinear  lift  range.  Derivation  of  the  unknowns  A, 

B,  C  and  N  is  described  in  Section  4.3.7. 

Two  other  user  options  are  available  from  the  routine;  (a)  the  user 

may  input  only  a  ,  or  (b)  the  user  inputs  only  a*.  Since  both  «  and  a* 
O  "  O  ” 

are  required  to  estimate  the  lift  variation  by  the  preceding  technique,  the 
subroutine  will  provide  an  estimate  for  the  missing  parameter  from  a  qua¬ 
dratic  expression.  Specifically,  a  quadratic  polynomial  can  be  faired 
through  the  nonlinear  lift  region  if  a  *  is  an  unknown.  Applying  the  gener¬ 
alized  boundary  conditions  to  a  polynomial  of  order  two,  and  solving  for  a* 
will  yield  an  estimate  for  this  unknown.  Conversely,  if  aQ  is  not  input,  it 
can  be  determined  in  a  similar  manner. 


■  Vf*  ’■'S'i 


I 


The  relationships  used  are  as  follows: 


1.  *  not  inp'Jt 

VV  +2l“o-“c, 


«o  not  input 


2  a  ^  a 

L  *  -  C. 

max  L 


CL  C 

max  L 


If  neither  a  nor  a*  are  user  inputs,  no  solution  is  possible,  but 
o  n 

the  program  calculated  values  for  and  ac  are 

a  max  L 

available  as  partial  output.  max 

4.3.2  Transonic  Wing  Drag  due  to  Lift,  Cn 

_  j 

The  programmed  procedure  for  computing  the  ratio  Cp  /C^  is  exactly  as 
described  in  Datcom  section  4. 1.5. 2.  The  method  does  a  three  dimensional 
table  lookup  for  Figure  4.1.5.2-55a  (A  tan  (A^g)  *  0)  and  for  Figure 
4.  1.5.  2-55b  (A  tan(ALg)  *  3).  Figure  4.1.5.2-55c  shows  a  linear  relation¬ 
ship  of  the  dependent  variable  (t/c)”*/^  Cq^/C^  as  a  function  of  the  tran¬ 
sonic  similarity  parameter  A  tan  (Ale)  f°r  each  value  of  the  ratio  (M^  - 
l)/(t/C)2/3;  it  was  assumed  that  this  linear  relationship  would  hold  for 
all  other  taper  ratios  other  than  0.50.  Therefore,  linear  extrapolations  on 
all  varibles  would  be  performed  if  required. 

This  method  was  programmed  in  subroutine  W1NGCL  with  the  calculation 
for  wing  Cl«  Since  Cl  is  required  to  calculate  C^,  the  calculation 
of  wing  Cl  would  enable  the  calculation  of  this  parameter  if  Cl  is  not 
input  as  experimental  data.  The  routine  will  not  overwrite  experimental  data 
input,  and  thus  the  user  oriented,  features  are  retained. 

The  r~tio  Cdl/Cl^  is  available  from  the  routine  and  will  be  output 
for  user  reference  if  Cq^  cannot  be  calculated. 

4.3.3  Transonic  Wing  Roll  Derivative,  C„ 

B 

Like  the  wing  Cj)L  calculation  described,  the  method  of  Datcom  Section 


5. 1.2.1  requires  wing  lift  to  calculate 


from  the  relationship  C^  /Cl, 

B  6 


equation  5.  1.2. 1— c.  Thus,  this  method  is  also  programmed  in  subroutine 
WINGCL.  The  calculated  value  f  j:  C,  will  not  overwrite  any  experimental 


IS--V  .**'.▼*.,'*  -t* 


data  input.  The  ratio  C^  /Cl  is  provided  if  the  calculation  for  C^  cannot 
be  completed.  No  exceptions  are  taken  for  the  Datcom  method.  The  ratio, 
C i  / C i  at  Mach  numbers  0.6  and  1.4  are  obtained  by  calling  the  subsonic 
and  supersonic  aerodynamic  modules. 

4.3.4  Transonic  Wing-Body  Roll  Derivative,  C„ 

The  derivative  C.,  will  be  calculated  by  Datcom  equation  5.2.2.1-d  if 

■  ,  *  6 

the  wing-body  lift  coefficient  variation  with  angle  of  attack  is  supplied,  or  • 

computed  as  described  above.  The  ratio  C^/Cl  is  given  as  partial  output 

if  the  lift  variation  is  not  specified.  This  method  is  implemented  exactly 

as  described  in  Datcom  and  is  programmed  in  subroutine  WBCLB.  Since  / 

,  <  ,  B 

Cl  at  Mfb  and  Mach  1.4  are  required  input  items  for  this  method,  they  are 

calculated  by  calling  the  appropriate  aerodynamic  modules. 

4.3.3  Transonic  Wing-Body-Tail  Drag  Coefficient,  Cj) 

This  method  is  a  "method  for  all  speeds"  as  described  in  Datcom  Section 

4.  3. 3. 2,  and  is  incorporated  in  exactly  the  same  manner  as  presently  pro¬ 
grammed  for  the  subsonic  solution.  This  method,  as  programmed  in  subroutine 
CDWBT,  require  the  following  experimental  data  inputs: 

1.  cDyg  vs  an8^e  attack 

2.  C[)H  vs  angle  of  attack 

3.  CLh  vs  angle  of  attack 

4.  q/q^  vs  angle  of  attack 

5.  «  vs  angle  of  attack 

6.  Cn  or  Cn 

°V  °WBT 

If  Cqo  is  not  an  experimental  data  input  item,  the  program  will 

calculate  it  from  the  estimated  Cn  calculated  as  follows: 

°WBT 

Cdov  *  Cd°wbt  Cd°wb  C°°h 

No  partial  output  is  available  from  this  method. 

4.3.6  Transonic  Wing-Body-Tail  Zero  Lift  Drag  Coefficient,  Cj>o 

This  method  follows  exactly  the  method  of  Datcom  section  4.5.3. 1, 
and  is  programmed  as  subroutine  WBTCDO.'  This  routine  does  not  require 
experimental  data  input,  although  experimental  data  input  is  an  optional 
feature  for  this  routine. 
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Utilizing  appropiate  configuration  description  parameters  the  prog  sm 
computes  the  drag  divergence  Mach  number,  Mu,  from  Figure  4.5.3,1-19.  The 
experimental  data  input  allows  the  user,  at  his  option,  to  select  -he  type  of 
general  configuration  to  be  used  in  computing  Mu.  The  three  options 


are .: 

o  A  -  Straight  wing  designs  without  area  rule, 
o  B  -  Swept  wing  designs  without  area  rule. 

o  C  -  Swept  wing  designs  incorporating  transonic  area  rule  theory. 

The  program  default  options  are  as  follows: 
o  No  wing  sweep  —  General  Configuration  A 

o  Swept  wing,  configuration  type  not  defined  -  General  Configuration  B 
The  general  configuration  types  are  defined  by  the  parameter  ITYPE, 
where  ITYPE=1  for.  configuration  type  A,  ITYPE=2  for  configuration  type 
B,  and  1TYPE=3  for  type  C.  In  the  case  of  configuration  type  C,  the  line 
for  type  C,  in  Figute  4.5.3.1-19,  was  linearly  extrapolated  and  programmed. 
All  extrapolations  in  this  figure,  with  the  exception  of  thickness  ratio,  are 
assumed  to  be  linear;  thickness  ratio  is  extrapolated  in  a  quadratic  fashion. 


With  Md  calculated  from  Figure  4.5.3.1-19, 

r  curve  across  the  transonic  Mach  regime. 
Do 

used  to  fair  the  curve: 


it  is  necessary  to  fair  the 
The  following  criteria  was 


A  polynomial  fairing  of  the  same  type  as  used  for  the  wing  nonlinear 
lift  coefficient  is  used  here  and  has  shown  acceptable  results.  , 

The  values  of  Cqo  at  Mach  .7  and  1.1  for  this  method  are  obtained  by 

calling  the  subsonic  and  supersonic  aerodynamic  modules. 
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4.  3.  7  Data  Fairing  Technique 

The  data  fairing  technique  used  for  computing  the  nonlinear  lift  region 
of  transonic  wings  and  the  transonic  wing-body-tail  zero  lift  drag  coef¬ 
ficient  was  chosen  for  its  powerful  features  and  ease  of  application. 
The  general  fairing  formula  is  a  polynomial  whose  form  is: 
y  =  A  +  B(X-Xq)  +  C(X-X0)N 

where  A,  B,  C  and  N  are  unknowns.  Given  the  values  of  y  and  dy/dx  at  two 
points,  Xy  and  Xj ,  four  simultaneous  equations  can  be  written.  These 
equations  solved  simultaneously  for  the  four  unknowns  yield  the  following 
results: 


A  =  y_ 


B  =  j*  0  X=X' 
dx  0 


yi  *  yo  (dx\  <xrV 

.  (XI  -  vN 

N.  [<&x,  -  (&x0]<V*o 


yi  "  yo  (dx;y  (xi  "  V 
o 


The  general  equation  reduces  to 


y  -  yo  +  €>XC  «-V  + 


[yI  *  yo  - 


(x 


'M  fiy 


This  equation  is  valid  for  Xq  <  X  £  Xj  and  (dy/dx)x0  ¥  (dy/dx)xj. 
Neither  of  these  conditions  is  violated  in  this  application.  The  range  of 
values  of  X  will  always  fall  between  Xq  and  X|  because  of  the  program 
logic,  and  in  the  nonlinear  lift  region  the  slopes  at  Xq  and  Xj  will 
never  be  equal.  For  the  transonic  wing-body-tail  versus  Mach 
fairing  the  Datecom  relation  (dC^/dM)  ■  0.10  at  M*Mp. 
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4.4  SUBSONIC  WING  Cm,  SUBSONIC  AND  SUPERSONIC  WING  AERODYNAMIC  CENTER,’  SUb- 
SONIC  WING-BODY  Cm,and  SUBSONIC  WING-BODY-TAIL  Cm 

The  subsonic  wing  pitching  moment  variation  with  angle  of  attack 
follows  Datcom  Method  1  of  Section  4. 1,4.3,  and  is  programmed  in  subroutine 
CMALPH.  The  method  is  applicable  to  those  configurations  whose  wing  aspect 
ratio  satisfies  the  following  criteria : 


a  <  6  ("LOW  ASPECT  RATIO") 

(l+c,)  cos  a  LE 

For  "high  aspect  ratio"  configurations,  the  default  wing  aerodynamic 
center  is  either  the  quarter-chord  of  the  wing  mean  aerodynamic  chord, 
or  the  user  input  value  (variable  name  in  the  planform  section  charac¬ 
teristics  namelists).  This  value  is  used  in  computing  pitching  moment  for 
the  wing  up  to  the  angle  of  attack  where  the  wing  lift  deviates  by  more  than 
7.52  from  the  linear  value;  at  this  point  the  method  is  no  longer  valid. 

There  are  no  methods  in  Datcom  or  Digital  Datcom  for  supersonic  wing 
pitching  moment,  though  the  wing'X^c  is  estimated  to  be  at  the  wing  plan-^ 
form  centroid  for  unswept  leading  edges,  and  computed  using  the  method  and 
design  charts  of  Datcom  section  4. 1.4.2  for  other  surfaces.  These  supersonic 
data  ate  computed  in  subroutine  SUPLNG. 

There  is  no  Datcom  method  for  computing  the  wing-body  pitching  moment  in 
any  Mach  regime.  Digital  Datcom,  however,  computes  the  subsonic  wing-body 
pitching  moment  using  the  following  formulation  (programmed  in  subroutines 
WBCMO  and  WBCM): 

o  Compute  (C0q}ub  from  regression  formulation  of  Datcom  Section 
4. 3.2.1,  programmed  in  WBCMO.  If  the  method  is  not  applicable, 
(CmQ)wB  computed  from  Method  1. 

o  Compute  the  wing-body  aerodynamic  center  from  Datcom  Section  4.3.2. 2 
(WBCM),  where  Equation  4.3.2.2-a  is  used  at  all  speeds. 

o  The  wing-body  Ca  curve  is  then  computed  as 

Stt  “  Sjwb  +  °»CL  +  C®CD 
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where  Cm  ,  is  the  pitching  moment  due  to  lift  obtained  by  Integra- 

la 

ting  the  curve  of  versus  from  *  0  and  to  at  the  desired 
angle  of  attack,  and  Cmc  is  the  pitching  moment  due  to  wing-body  drag 
located  at 

Subsonic  wing-body-tail  pitching  moment  versu_  angle  of  attack  is 
computed  by  Digital  Datcom  in  subroutine  WBTAIL,  though  there  is  no  Datcom 
method  for  this  parameter.  The  method  formulation  used  is  as  follows: 


j  H  j  WBT  j  WB 


(SL-(%)vb+ 


(x  -  x  ).. 

ac  eg  II 


C  )  COS  (a) 

/  U  J 


(cDj)h  (Vfl.)j  SIN  (O).  +  (.V..-it-  (cDj)b  (1/<W)j  COS 


i\  )„ 


SIN  (a) 


4.5  TRANSONIC  BODY  CL  FAIRING  AND  TRANSONIC  BODY  C^  FAIRING 
•  _ ; . . .  a _ _ 

The  transonic  C^  and  CD^  derivatives  for  the  body  alone  configura¬ 
tion  is  interpolated  linearly  between  thu  subsonic  (M  *  0.60)  and  supersonic 
(M  •  1.40)  Mach  regimes  in  subroutine  B0DVRT. 
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4.6  SUBSONIC  ASYMMETRICAL  BODY  CL,  SUBSONIC  ASYMMERICAL  BODY 

Cm,  AND  SUBSONIC  ASYMMETRICAL  BODY  CDf>,  Cp _  ' 

Digital  Datcom  body  solutions  generally  include  lift,  drag,  and  pitching 
momenc  coefficients.  In  the  transonic  speed  regime  the  solutions  are  re¬ 
stricted  to  lift  and  pitching  moment  slopes,  and  drag  coefficients. 

4.6.1  Subsonic  Bodies 

Subsonic  body  analysis  computes  lift,  drag,  and  pitching  moment  coef¬ 
ficients  for  either  axisymmetric  or  cambered  bodies.  Digital  Datcom  body 
methods  are  identical  to  Datcom  except  for  the  base  drag.  Digital  Datcom 
calculates  base  drag  using  a  minimum  base  area  equal  to  30%  of  thq  body 
maximum  cross-sectional  area. 

The  cambered  body  pitching  moment  method  is  not  defined  in  Datcom 
and  is  therefore  described  in  detail.  For  clarity,  the  lift  method,  which  is 
defined  in  Datcom,  is  also  described.  These  body  methods  (subroutine 
B0D0PT)  are  executed  when  the  parameters  Z\j  and  Zg  are  user  specified 
(namelist  B0DY).  The  method  predicts  the  zero  lift  angle  of  attack, 
zero  lift  pitching  moment,  and  body  lift  and  pitching  moment  versus  angle  of 
attack.  The  Datcom  drag  methods  are  retained. 

Zero  lift  angle  of  attack  and  pitching  moment  are  calculated  utilizing 
conventional  mean  line  theory.  The  equations  are: 
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These  parameters  are  defined  in  Figure  25. 

Lift  and  moment  for  asymmetric  bodies  are  calculated  by  employing  a 
modified  version  of  Polhamus's  leading-edge  suction  analogy  (References 
2  and  3).  Polhamus  considers  two  components  of  lift,  a  potential  flow 
term,  C^,  and  a  vortex-lift  term  Clv.  Both  of  these  terms  are  a 
function  of  body  aspect  ratio  (A)  and  are  defined  as  follows: 

CL  “  cLp  +  CLv 

CLp  *  Kp  sin  a  cos'-  a 
*  Ky  sin^  a  cos  a 
a  *  angle  of  attack 
Kp  and  K y  are  obtained  from  Figure  26. 

The  Polhamus  vortex  lift  equation  i^ust  be  modified  to  make  it  applicable 
to  thick  bodies  because  the  onset  of  vortex  lift  for  such  configurations  is 
not  at  zero  angle  of  attack  as  it  is  with  flat  plate  wings.  The  thick  body 
angle  of  attack  for  onset  of  vortex  lift  (av)  can  be  correlated  with  the 
fineness  ratio  (FR)  >.nd  tne  thickness  ratio  (TR.)  of  the  body  as  shown  in 
Figure  27a.  The  body  thickness  parameters  are  shown  in  Figure  27b.  Experi¬ 
mental  data  used  in  correlation  are  presented  in  References  4  through  7.  The 
redefined  lift  expressions  for  thick  bodies  are  as  follows: 

Clp  ■  Kp  sin  a  cos2  u 

C'Ly  ■  Ky  sin2  (  a  -  oi  v)  Cos  ('  a  -  a’y) 

c'l-cLp  +  c*Lv 

The  body  pitching  moment  is  obtained  by  estimating  the  center-of- 
pressure  locations  of  both  the  potential  and  vortex  lift  components. 
The  total  pitching  moment  is  equal  to  the  sum  of  the  moments  produced 
by  the  lift, forces  acting  at  their  respective  center-of-pressure  loca¬ 
tions  plus  the  zero  lift  pitching  moment.  The  potential  lift  center-of- 
pressure  location  employed  stems  from  slender  body  theory  and  is  presented  in 
Figure  28  as  a  function  of  n.  The  equation  for  the  powerlaw  planform  is  of 
the  form  R  *  Rmax  (X/L)n.  The  program  computes  an  exponent  n  that  closely 
approximates  the  input  planform  area.  The  potential  lift  center-of-pressure 
location  Is  obtained  from  Figure  28  or  the  equation, 

Xcp/L  -  2n/(2n+l> 
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FIGURE  25  ASYMMETRIC  BODY  GEOMETRY  INPUTS 


DEGREES 


FINENESS  RATIO  -  FR 
OATCQM  FIGURE  4.2.1.2-37 

FIGURE  27a  CORRELATION  OF  ay 


BODY  PLANFORW 


FIGURE  27b  BODY  THICKNESS  PARAMETERS 
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of  pressure  Is  assumed  to  be  located  at  the  total  planform 
The  equation  for  the  body  pitching  moment  coefficient  is: 

^ra  *  ^mQ  +  ^mp  +  cmy 
Cmp  *  <*CG  -  XCP)/^ 

C*v  “  CNV  <*CG  “  X)/L 

where  X  is  the  location  of  the  total  planform  center  of  area  measured  from 
the  body  nose.  The  method. is  applicable  at  angles  of  attack  equal  to  or 
greater  than  the  wing  maximum  lift  angle  of  attack.  , 

4.6.2  Transonic  Bodies 

Digital  Uatcom  body  solutions  are  restricted  to  lift  and  pitching 
moment  slopes,  and  drag  coefficients  in  the  transonic  speed  regime.  These 
data  are  computed  by  performing  a  linear  interpolation  between  the  subsonic 
(M  ■  0.60)  and  supersonic  (M  **  1.4)  Mach  regimes. 

Subroutines  that  implement  the  transonic  body  methods  are  B0DYRT, 

SUPB0D,  TRS0NI,  and  TRS0NJ. 

4.6.3  Supersonic  Bodies 

Supersonic  body  analysis  provides  solutions  for  lift,  drag  and  pitching 
moment  coefficients.  Datcom  methods  for  lift,  pitching  moment  slope, 
and  drag  coefficient  require  the  body  to  be  synthesized  from  a  combina¬ 
tion  of  body  components  comprised  of  a  nose,  after-body,  and/or  tail  seg¬ 
ments.  Digital  Datcom  requires  synthesized  body  configurations  to  be  either 
nose  alone,  nose-after  body,  nose-after  body-tail,  or  nose-tail  segment 
combinations. 

Some  of  the  Datcom  body  drag  methods  in  this  speed  regime  have  not  been 
implemented  in  Digital  Datcom.  The  effects  of  blunted  noses  on  drag  are  not 
incorporated  since. the  body  lift  and  pitching  moment  methods  do  not  reflect 
the  Influences  of  this  parameter.  Some  of  the  Datcom  interference  drag 
methods  are  also  deleted.  In  this  case,  the  methods  were  omitted  because  of 
their  limited  range  of  applicability. 

Calculation  of  wing-body,  or  horizontal  tail-body,  stability  requires 
the  lift  curve  slope  of  the  body  ahead  of  the  wing  or  horizontal  tail.  Body  . 
Ctf  methods  are  executed  for  the  portion  of  the  body  ahead  of  the  wing,  if 
the  wing  is  present;  the  portion  of  the  body  ahead  of  the  horizontal  tail,  if 
the  horizontal  tail  is  present;  and  entire  body. 


Vortex  lift  center 
centroid  of  area. 
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All  methods  are  implemented  by  subroutine  SUPB0D  except  for  a  portion 
of  the  drag  methods  contained  in  subroutine  F1G26. 

4.b.4  Hypersonic  Bodies 

Hypersonic  body  analysis  is  performed  at  user  designated  Mach  numbers 
that  are  equal  or  greater  than  1.4.  Iri  this  speed  regime,  Digital  Datcom 
stability  solutions  include  lift,  drag  and  pitching  moment  coefficients. 

Hypersonic  body  analyses  for  lift  and  pitching  moment  slopes  and  drag 
coefficients,  like  the  supersonic  body  methods,  require  the  body  to  be 
synthesized  from  a  combination  of  body  segments.  Hypersonic  body  analysis  is 
unlike  the  other  Datcom  hypersonic  configuration  analyses  since  the  methods 
are  defined  independent  of  the  supersonic  results.  Body  for  portions 
of  the  body  ahead  of  the  wing  and/or  horizontal  tail  are  also  calculated. 

The  methods  are  implemented  in  subroutine  HYPBOD.  A  small  portion 
of  the  drag  methoas  are  found  in  subroutine  FIG26. 

4.7  TRANSONIC  WING-BODY  CL 

The  transonic  wing-body  lift  coefficient,  if  not  input  using  name- 
list  EXPR — ,  is  computed  in  subroutine  WBCLB  using  the  following  equations: 


In  computing  the  transonic  wing-body  pitching  mot  .nt  slope,  the  center 
of  pressure  of  body-wing  carryover  is  linearly  interpolated  between  the 
values  obtained  at  Mach  0.60  and  Mach  1.40  in  subroutine  TRANCM. 

4. 8  WING-BODY-TAIL  MOVEABLE  HORIZONTAL  TAIL  TRIM 

Tne  all  moveable  horizontal  tail  incidence  required  to  trim  the  vehicle 
(Cm  ,  ■  0)  at  angle  of  attack  is  calculated  in  subroutine  TRIMR2.  At 

trim,  the  forces  on  the  tail  are  and  (trimmed  lift  and  drag),, 
and  are  referenced  to  the  local  flow  at  a  tail  angle  of  attack  of  (  a  -  <H). 
Since  these  trimmed  forces  are  located  at  the  tail  aerodynamic  center,  which 
is  known,  the  total  body  moments  can  be  summed  as  follows ; 


can  be  expressed  as 


Hence,  the  only  unknown  is  ,  the  tail  lift  at  trim,  which  can  be 
evaluated.  From  Sketch  (a)  note  that 


V 


£cn. 


’  trim 


VIEW  IN  PLANE  OF  SYMMETRY 
a  =  Airplane  ar.glc  of  attack  (positive  as  shown) 

xH  =  Distance  from  c.g.  to  quarter-chord  point  of  horizontal- stabilizer  MAC 

fl  =  Angle  defined  by  intersection  of  xH  with  FRP  (positive  as  shown 
with  horizontal  stabilizer  above  c.g.) 


Sketch  (a) 
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The  moment  equation  reduces  to 


qu  qu  . 

cm  +  cm  “  -  C.  ~  —  cos  (n  -  a  +  Em) 

"wb  "oh  q-  lh  q.  H 


+  Cn  *  -ip— -  —  —  sin  (fi  -  a  +  eH)  3  0 

°0H  'AHeH  q-  Ct,  H 


Letting  6  =  fi  -  a  +  £„  and  re-arranging  yields  a  quadratic  on  C 


qH  XH  ,  &  (CLH) 

^  r  5  »vH 

w 

.  Ji  cos  6  (c,  ) 

q°°  cw  H 


qH  XH  • 


^D0H  q~  c,, 


sin  5  +  C„  +  C 


^WB  M0H  q° 


Simplifying, 

qH 

^  (C,  )2  -  C,  *  Cn  tan  a  t  ^ 

ffAHeH  LH  LH  °0H  %  cos 

9®  Tf 
W 

From  the  quadratic  formula. 


In  this  form,  the  equation  becomes  invalid  for  *  0,  and  can  be  further 
reduced  to 


Qu 

r  +  r 

M0H  q°° 

XU  % 

~  ~  COS  6 

L  cw  “ 


+  C_  tan  6 
d0H 


1 — 

- 

r 

h 

tan  6 

,AHeH 

.  '  J 

r  +  r 

“wb  mqh  %.  +r 
XH  qH  LD, 

—  -2-  COS  8 

c  q» 

uw 


tan  5 


OH 


A  plus  sign  in  front  of  the  radical  is  the  valid  solution,  otherwise  at 
6*0  the  solution  is  undefined.  This  result  is  similar  to  Datcom  equation 
4.5. 3. 2-e,  with  the  exception  of  the  term 

Once  the  tail  lift  at  trim  (CLh)  has  been  determined,  a  variation 
of  Datcom  equation  4.5.  1.2-a  can  be  used  to  calculate  the  tail  incidence 


iH* 


CL„  =  ClH  (KH(B)  +  KB(H)) 


-H 


+  C*  (a.)  [kH(B)  +  kB(Hj] 
aH  H 


+  I, 


'B(K) 


'  '  H 


(b/2  -  b*/2) 

TO 


CL*  °eff 
aH 


where  i3  the  pseudo  lif t-curve-slope  of  the  horizontal  tail  in  the 

presence  of  the  body, 


C,* 


L*  =  CL  <KH(B)  +  XB(H)) 
°H  °H 


103 


Cl^'  and  C^'  are  the  horizontal  tail  lift  and  lift  curve  slope  at 

(a  -  cH  +  aQH) 

and  aaff  is  the  effective  angle  of  attack  of  the  horizontal  tail  in  the 

presence  of  the  body  /.  .  \ 

I  *H(I)  +  kB(H)  \ 

“eff  *  «H  *  “OH  % 


The  incidence  angle  to  trim  can  then  be  solved  directly,  and  becomes 


£H  +  °0H  )  !Vr 


C  S)  '  *  "sw)  [V  +  CCh  (°  ‘  tH  +  “0H)  (svf)  ( 2|>«  }\ 

H(B)j  L  aH  VB(H)  \2itVr  /  _H  y  b/2  /  LotH 


Once  the  tail  lift  and  drag  at  trim  has  been  computed  the  panel  hinge 
moment  about  the  pivot  point  can  also  be  computed.  Since  Cj^j  and  Cj)H  are 
are  referenced  to  the  local  flow,  they  must  be  computed  relative  to  the 
freestream  flow.  Relative  to  V^,  ,  trim  lift  and  drag  are 

Ru 

C.  =  (C.  COS  e  ■?  Cn  Sin  e)  — 

Vim  ht  ht  q- 


TrAHeH 


The  pitching  moment  trimmed  is 


-Hh  Lh  c 

"TRIM  "TRIM 


COS  6  +  C 


sin  6 


The  hinge  moment  about  the  pivot  point  is 


CHM  3  CLU  C0S  a  +  CDU  Sin  a 

_  "trim  "trim 
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4.9  WING-BODY-TAIL  TRIM  WITH  CONTROL  DEVICES 


Conf iguratior.  trim  with  wing  or  horizontal  tail  control  devices  is 
performed  in  subroutine  TRIMRT.  The  method  programmed,  which  is  not  a 
Datcom  method,  essentially  does  a  table  look-up  of  the  control  device 
incremental  pitching  moment  coefficient  versus  control  deflection  for 
the  deflection  required  to  trim.  The  incremental  lift  coefficient  and 
drag  coefficient  ire  then  obtained  by  performing  table  look-ups  for  these 
variables  (which  are  a  function  of  control  deflection  angle)  at  the  trimmed 
control  deflection. 

4. 10  STANDARD  ATMOSPHERE  MODEL 

Incorporation  of  a  standard  atmosphere  model  (subroutine  ATMOS)  into 
Digital  Datcom  provides  input  and  output  flexibility  for  the  user.  The 
program  can  operate  on  Mach  number  and  altitude  as  separate  independent' 
variables.  The  addition  of  vehicle  weight  and  flight  path  angle  permit 
calculation  of  equilibrium  flight  conditions. 

The  program  allows  the  user  to  input  either  Mach  number  or  velocity 
as  an  independent  variable  for  speed  reference.  If  velocity  is  input, 
the  free  stream  static  temperature  must  be  available  so  that  Mach  number 
can  be  calculated.  The  user  will  also  have  the  option  to  specify  a  flight 
altitude,  or  static  pressure  and  temperature,  as  an  independent  variable 
defining  the  atmospheric  conditions.  If  altitude  is  specified,  pressure  and 
temperature  will  be  found  using  the  "U.S.  Standard  Atmosphere,  1962." 

The  user  may  input  up  to  20  Mach  or  velocity  points.  If  Mach  number  is 
input,  the  veloqity  will  be  calculated  for  each  point  where  atmospheric  data 
are  input.  When  velocity  is  input  the  Mach  number  will  be  calculated  using 
atmospheric  conditions.  If  velocity  is  input  instead  of  Mach  numbers  and 
atmospheric  conditions  are  not  defined,  an  error  message  will  be  written  and 
Mach  numbers  will  be  calculated  using  a  speed  of  sound  of  1000  ft/sec. 

The  user  may  also  input  up  to  20  atmospheric  conditions.  The  atmosphere 
may  be  defined  by  altitude,  pressure  and  temperature,  or  Reynolds  number.  If 
the  altitude  is  given,  pressure  and  temperature  will  be  determined  using  the 
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atmosphere  model  developed  in  Reference  9.  The  Reynolds  number  will  be 
calculated  using  the  following  equation  (in  the  foot-pound-second  system  of 
units) : 

RN/L  -  pV/y  -  1.2527  X  106  PM  (T  +  198.6)/T2 
This  equation  was  derived  using  the  following  relationships: 
p  -  P/RT 

v  *  m  Art  ' 

P  -  2.270  X  10~8  T1 *5/lT  +  198.6) 

If  the  Reynolds  number  is  not  input  and  cannot  be  calculated,  an  error  mes¬ 
sage  will  be  written  and  the  Reynolds  number  will  be  set  to  5  x  10^/ft. 

Given  the  vehicle  weight,  flight  path  angle,  and  atmospheric  condi¬ 
tions,  the  equilibrium  flight  aerodynamic  data  can  be  determined.  Equilib¬ 
rium  flight  is  achieved  when  the  following  relationship  is  satisfied. 

WT  *  (CL  cos  6  -  Cp  sind  )  qS 

Along  with  the  untrimmed  aerodynamic  output,  the  level  flight  (6*0) 
lift  coefficient  will  be  output.  Trim  data  output  will  provide  an  addi¬ 
tional  line  of  output  at  the  equilibrium  flight  conditions  (subroutine 
FLTCL). 
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SECTION  5 


SYSTEM  RESOURCE  REQUIREMENTS 


Digital  Datcom  is  a  large  and  rather  complex  computer  program  which 
requires  specific  computer  resources  to  execute  within  a  fixed  core  require¬ 
ment.  The  program  is  written  to  conform  to  the  American  National  Standards 
Institute  (ANSI)  Standard  Fortran  IV.  Certain  computer  resources  must  be 
available  to  make  the  program  operational  without  modifications  These 

resources  are: 

o  Six  disk  files  or  scratch  tapes  are  required  for  manipulation 
and  retrieval  of  input  data.  The  logical  I/O  units  used  are  8,  9, 
10,  11,  12  and  13.  These  logical  units  are  in  addition  to  logical 
unit  5  (read)  and  unit  6  (write). 

o  The  system  must  have  capability  for  primary  and  secondary  overlay 
structures. 

o  The  system  must  have  a  Fortran  compiler  which  provides  for  NAMELIST 
input  and  output,  and  statement  transfer  when  an  end  of  file  is 
detected. 

Each  logical  unit  referenced  by  the  program  is  reserved  for  a  specific 
purpose.  The  units  referenced  and  their  use  in  the  program  are  listed 
below: 


Unit  Program  Usage 

5  Standard  system  input  (card  reader) 

6  Standard  system  output  (printer) 

8  Storage  of  experimental  data  namelists  for  the  case  being 

executed 

9  Storage  of  input  namelists,  dxcept  experimental  data,  for  the 
case  being  executed 

10  Storage  of  experimental  data  namelists  for  a  single  Mach 

number 

11  Storage  of  all  input  data  after  processing  by  the  input  diagnos¬ 
tic  analysis  module  (C0NERR) 

12  Storage  of  extrapolation  messages  for  processing  by  overlay  57 

13  Storage  of  output  data  for  use  with  the  Plot  Module  as  a  post¬ 
processing  option 


SECTION  6 


PROGRAM  CONVERSION  MODIFICATIONS 


6. 1  GENERAL  REMARKS 

The  program  was  written  in  Fortran  IV  for  the  CDC  Cyber  175  computer 
system.  Several  program  modifications  may  be  required  to  run  under  other 
Fortran  compilers  or  computer  systems.  It  is  recommended  that  users  imple¬ 
menting  the  program  for  their  computer  system  become  familiar  with  their 
installation  operating  system  and  Fortran  compiler  requirements.  Users  are 
forewarned  that  program  core  requirements  and  run  times  discussed  in  this 
report  may  no  longer  be  valid. 

6.2  PROGRAM  STRUCTURE 

The  program  is  composed  of  a  root  segment  overlay  (overlay  0),  f iffy- 
seven  primary  overlays  and  twenty-eight  secondary  overlays.  Table  7  shows 
the  overall  program  structure  and  lists  those  routines  that  are  contained  in 
each  overlay.  In  the  CDC  system,  the  first  routine  in  an  overlay  is  called  a 
"program"  and  subsequent  routines  "subroutines.”  Several  subroutines  appear 
in  more  than  one  overlay.  These  subroutines  are  called  "common  decks"  and 
are  listed  in  Table  8. 

6.2.1  Calls  to  Overlay 

All  primary  overlays  are  called  by  the  root  segment  overlay,  and 
secondary  overlays  called  by  their  respective  primary  overlay  using  the 
calling  sequence 

CALL  OVERLAY  (4LDATC,  XX,  YY,  6HRECALL) 

where:  DATC  is  the  disc  file  where  the  overlay  is  located,  . 

XX  is  the  primary  overlay  number  in  decimal,  and 
YY  is  the  secondary  overlay  number  in  decimal. 

Hence,  each  overlay  is  written  to  a  disk  file  with  the  name  “DATC." 

Users  should  refer  to  the  Fortran  reference  manual  for  their  system  and 
determine  the  correct  overlay  calling  procedure. 

6.2.2  Common  Decks 

Several  subroutines  are  used  in  more  than  one  overlay.  The  most  commonly 
used  routines  are  located  in  the  root  segment  for  access  by  all  overlay 
programs  and  subroutines.  However,  several  decks  are  used  by  only  a  few 
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routines  and  placing  them  in  the  root  segment  would  require  an  increase  in 
overall  program  core  size.  In  order  to  maintain  a  low  core  requirement, 
these  common  decks  are  located  in  each  overlay  in  which  it  is  referenced. 

Warning  -  Not  all  systems  allow  two  routines  to  have  the  same  name  even 
though  they  are  identical.  If  the  user's  system  does  not  allow  this  option, 
three  alternatives  are  available  as  follows: 

o  Rename  each  deck  that  is  common,  and  change  the  calling  sequence  to 
it. 

natives  are  available  as  follows: 
o  Place  all  common  decks  in  the  root  segment  (overlay  0)  and  re¬ 
move  the  deck  from  each  associated  overlay.  The  user  will  increase 
the  overall  program  core  requirement  by  using  this  technique,  how¬ 
ever,  it  is  easier  than  the  procedure  outlined  above, 
o  On  some  systems  that  have  multiple  region  capability,  these  common 
decks  can  be  placed  in  a  separate  overlay  region. 

6. 2. 3  "OVERLAY"  and  "PROGRAM"  Cards 

Each  primary  and  secondary  overlay  main  program  contains  these  two 
cards.  The  CDC  Fortran  compiler  requires  all  overlays  to  begin  with  an 
“OVERLAY"  card  followed  by  a  main  program  which  begins  with  a  "PROGRAM"  card. 
These  must  be  replaced  by  corresponding  code  required  by  the  operating  system 
and  compiler  being  employed. 

6.2.4  End  of  File  Tests 

Routines  INPUT,  C0NERR  and  XPERNM  utilize  a  transfer  on  end  of  file. 

This  statement  must  be  modified  for  the  Fortran  compiler  being  used. 

6.2.3  Use  of  *' U.IJSED"  and  "KAND" 

These  constants  are  set  in  BL0CK  DATA.  The  value  for  "UNUSED"  is 
set  in  the  program  as  10  60.  It  is  sometimes  used  as  a  program  flag  and  is 
used  for  initializing  all  variable  arrays  to  some  number  other  than  zero. 

The  value  for  "UNUSED"  can  be  changed  if  desired  and  must  be  defined  in  BL0CK 
DATA  as  a  small  positive  number.  The  variable  "KAND"  defines  the  alphabetic 
character  used  by  the  NAMELIST  1’iputs.  It  is  set  to  '$'  for  CDC  systems. 
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SECTION  7 


PROGRAM  DECK  DESCRIPTION 

This  section  contains  a  description  of  all  routines  in  Digital  Datcom. 
Table  7  lists  the  decks  by  overlay,  Table  8  lists  those  "common  decks"  in  the 
program,  and  Table  9  describes  the  purpose  of  each  deck  and  the  overlays 
referenced.  For  convenience.  Table  9  lists  the  routines  in  alphabetical 
order.  Table  10  discusses  the  use  of  each  of  the  variables  In  the  Digital 
Datcom  control  data  olocks.  The  description  of  the  plot  module  routines  is 
provided  in  Volume  III  of  this  report. 

A  complete  program  listing,  which  includes  Digital  Datcom  and  the 
Pxot  Module,  is  provided  as  a  microfiche  supplement  to  this  report.. 
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TABLE  7  DIGITAL  DATCOM  OVERLAY  DESCRIPTION 
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TABLE  7  DIGITAL  DATCOM  OVERLAY  DESCRIPTION 
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Deck  Name 
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ALI 

7,  10,  20,  28,  35 

ANGLES 

2,  13,  15,  16,  18 

ARCC0S 

12,  21,  41,  42,  50,  53 

ARCS IN 

21,  41,  42,  44 

B0D0WG 

7,  10,  20,  28,  35 

CALCALC 

31,  33 

CALCAO 

15,  16 

CDRAG 

3,  5 

CLMXBS 

,  15,  16 

CLMXBl 

24  (Both  Secondary  Overlays) 

CMALPH 

31,  33 

DFLC0N 

41,  53 

DMPARY 

11,  39,  42, ,46,  47,  49 

EQSPCE 

4,  6 

EOSPC1 

4,  6 

FIG53A 

3,  5 

FIG68 

21,  42 

GETMAX  . 

4,  6,  29 

INFTGM 

2,  21 

LIFTCF 

15,  16 

PRSCID 

12,  39,  42,  46,  47 

SETUP 1 

2,  18 

SIMUL2 

38,  42,  47 

SWRITE 

12,  39 

TABLEC 

7,  20,  25 

TABLES 

7,  24 

TBSUB 

7  >  24 

TBSUP 

7 1  24 

TBTRN 

7,  24 

TEST 

1,  ,34 

TLIN4X 

17,  25,  26,  52 

TLIPIX 

43,  44,  45,  54 

TLIP2X 

43,  44,  45,  54 

TLIP3X 

43,  44,  45,  54 

TRANF 

24  (Both  Secondary  Overlays) 

TRAPZ 

4,  6,  9,  19,  23,  26,  29,  37, 

WBCDL 

7,  24 

WBCMO 

7,  20,  25 

WBCM1 

25  (Both  Secondary  Overlays) 

WTGE0M 

2,  18 

WTLIFT 

15,  16 

YUP 

43,  44,  45,  54 

ZERANG 

1,  2,  13,  18 
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